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We compute the Next-to-Next-to-Leading Order (NNLO) QCD corrections to the partonic re- 
action that dominates top-pair production at the Tevatron. This is the first ever NNLO calculation 
of an observable with more than two colored partons, and/or massive fermions, at hadron colliders. 
Augmenting our fixed order calculation with soft-gluon resummation through Next-to-Next-to- 
Leading Logarithmic (NNLL) accuracy, we observe that the predicted total inclusive cross-section 
exhibits a very small perturbative uncertainty, estimated at ±2.7%. We expect that once all sub- 
dominant partonic reactions are accounted for, and work in this direction is ongoing, the perturbative 
theoretical uncertainty for this observable could drop below ±2%. Our calculation demonstrates 
the power of our computational approach and proves it can be successfully applied to all processes 
at hadron colliders for which high-precision analyses are needed. 



INTRODUCTION 



Due to the number of special features it possesses, the 
top quark has become one of the pillars of the physics pro- 
grams at the Tevatron and LHC. First, it has a uniquely 
large coupling to the Higgs sector of the Standard Model 
(SM). Second, the top quark is often the preferred decay 
mode of new massive objects, like heavy resonances, pre- 
dicted in many models of new physics. Third, it decays 
before it hadronizes which makes it possible to study top 
quarks without encountering the non-perturbative effects 
that obscure the production of lighter quarks. Finally, 
and notably, one of the most significant current devia- 
tions from the SM is in the top quark forward-backward 
asymmetry ApB HI- 

Motivated by these observations, in this work we make 
the first step towards a comprehensive increase of the pre- 
dictive power of the Standard Model in the top sector. 
Specifically, we evaluate the Next-to-Next-to-Leading 
Order (NNLO) QCD corrections to the total inclusive 
top-pair cross-section for the reaction dominating at the 
Tevatron, qq ^ ti -\- X. Our Tevatron prediction, based 
on this first ever NNLO calculation for a hadron collider 
process involving more than two colored partons and/or 
massive fermions, is almost twice as precise as the cur- 
rently available experimental measurements. We believe 
that our results are a strong motivation for further ex- 
perimental improvements in top physics. 

We expect to extend our work with fully differential 
results for top-pair production at the Tevatron and the 
LHC, as well as production of dijets, M^-l-jet, Higgs-|-jet, 
etc. All these processes will be instrumental in the ongo- 
ing and future searches for new physics and for assessing 
the workings of the Standard Model at the finest level. 



TOP PRODUCTION AT HADRON COLLIDERS 

The total top-pair production cross-section reads: 
crtot=V/ d/3$y(/3,M^)CT,J(/3,m^^2)^ 

where i,j run over all possible initial sate partons, 
is the partonic flux which is a convolution of the densi- 
ties of partons i,j and includes a Jacobian factor. The 
dimensionless variable /3 — y'l — p, with p = Arn^/s, is 
the relative velocity of the final state top quarks having 
pole mass m and produced at partonic cm. energy yfs; 
fi stands for both the renormalization (pB.) a-nd factor- 
ization scales (/ii?). 

For fip — fJ-R — fJ- the partonic cross-section reads: 
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where L = ln (/i2/m2) and as is the MS coupling renor- 
malized with Nl — 5 active flavors at scale fi^. The 
functions a^J^'"^^^ depend only on /3. 

The dependence on /ijj ^ can be trivially restored 
in Eq. ([2]). The LO resuh reads (J^jq = 7r;3p(2 + p)/27. 
The NLO results are known exactly [2|. The scale con- 

f2 ('2 2^ 

trolling functions cr ' ' and alj ' ' can be easily computed 

from the NLO results o'lj'', and can be found in Q. 

In this work, for the first time, complete results for the 
function are derived. Work towards the calculation 
of the remaining reactions ij — gg,gq,qq' is underway. 

(2) (2) 

We recall that the only information about and agg 
available so far was from their /3 — > limits [4] . 
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PARTON LEVEL RESULTS 
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In this paper we calculate the NNLO correction a^q 
to the partonic reaction qq — > tt + X. The result reads: 
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(/3) = Fo(/3) + Fi(/3)iVi + F^im. 
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The dependence on the number of light flavors N]^ in 
Eq. ([3]) is explicit. The function F2 is derived exactly: 
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The functions F, = F,^^^ 
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i = 0, 1, read: 
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F^^'> = crj"^ [0.0228463//32 + (-0.0333905 + 0.342203^ 



af^ [(0.0116822 
+0.353207L^ - 
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0.0277778L/5)//3 

0.240169L|] , (5) FIG. 1: The functions 
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and 10 X f[^'^^ (rescaled for 
improved visibility) as: a) discrete sets of calculated values, 
with errors, on the grid of 80 points (red and blue error bars) 
and b) the analytical fits given in Eqns. H8l7p (black lines). 
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9.18183/3^)^ 
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(0.90756/3 - 6.75556/3^ 4 
+ (-0.99749/3 + 27.7454^2 _ i2.9055/33)p2 
+ (-0.0077383^- 4.49375/3^ + 3.86854/3^)^3 
-0.380386/3V* + Lp(1.3894p + 6.13693/9^ 
+8.78276p3 - 0.0504095/) + 0.137881p , 

(-2.32235/3 + 44.3522/3^ - 24.6747^^)^ 
+ (2.92101/3 + 224.311/3^ + 21.5307/3^)^^ 
+ (2.05531/3 + 945.506/3^ + 36.1101^^ 
-176.632/3'*)p3 + 7.68918^V + Lp(3.11129p 
+100.125/3^ + 563.1/9=^ + 568.023^"*) , 
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where Lp = \n.(p) and = ln(/3). The functions F-f^^ 
constitute the threshold approximation to crg|^ Q mul- 
tiplied by the full Born cross-section cr^^-* and with the 

constant C^^^ (as introduced in Q) set to zero. 

The functions i*i are computed numerically on a grid 

(fit) 



of 80 points in the interval /3 € (0, 1). The functions F^ 
are then derived as a fit to the difference Fi — F, 



On Fig. [T] we present the fitted functions F^^^^ together 
with the calculated values for Fi — F^'^' in all 80 points, 
including the estimated numerical errors for each evalu- 
ation point. We note that the precision of our result is 
not limited by the quality of the fit, but rather by the 
precision of the numerical evaluation of -Fi^o- The abso- 
lute error on (J^qq , for Nl = 5, is bounded by 3.6 x 10^^. 
At the Tevatron this translates into a relative error on 
the cross-section of 3 x 10"'', which is negligible. 

The fits become unreliable very close to the high- 
energy limit /3 — >■ 1, i.e. beyond the last computed point 
/380 = 0.999. While this loss of precision is completely 



immaterial for top production (/38o = 0.999 is approx- 
imately the value of Praax for the LHC at 8 TeV), it 
might be an impediment for the description of lighter 
quark production, like bottom. The reason is that for 
very light quarks the partonic flux $ becomes strongly 
peaked towards /? « 1 which makes the hadronic cross- 
section sensitive to the behavior of the fits in this limit. 

From Eqn. ^ we can extract an approximate value 
for the two-loop constant term C^|\ as defined in 0], 

which translates into the hard matching constant H, 
as defined in Ref. We get the following values: 
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1195.82 -44.1841/VL-4.28168iV£, (9) 
84.8139 (for iVi = 5) . (10) 

Following the findings of Refs. @, we caution about the 
accuracy of the extraction of C^^^ (and therefore ) . 
Assuming a polynomial in /3 form for the fits F\^ ^ , we 
can extract C^q with a precision better than 10% (which 

implies £ (80,90)). This uncertainty has sub-permil 
numerical effect for top production at the Tevatron. On 
the other side, we note that if we allow into the fits terms 
containing powers of ln(/3), then C^q cannot be extracted 
with any reasonable precision (the reason being the finite 
size of the grid) . At any rate, the overall smooth behavior 
of the fits suggests that our extraction is reliable. 

It is interesting to compare the exact partonic cross- 
section Eq. ([3]) with the approximately known one Q. 
To that end on Fig. [5] we plot the partonic cross-section 
multiplied by the partonic fiux for Tevatron (i.e. the in- 
tegrand of Eq. ([TJ) for the following three cases: a) exact 
NNLO ©, b) approximate NNLO defined by setting F2 
and fI^^^ in Q to zero and c) as in b) with the additional 
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Approx NNLO, Exact Bom 
Exact NNLO 




partonic NNLO x Flux 
qq -> tt+X 

MSTW2008NNLO(68c. I. ) 



FIG. 2: Partonic cross-section times the partonic flux for 
Tevatron (same as the integrand of Eq. for three different 
partonic cross-section: a) exact NNLO (black thick line), b) 
NNLOapprox with exact Born aj^^ (blue dashed line) and c) 
NNLOa 



with leading Born a^^' 1/3^0 (red thin line) 



replacement a^j^ 



(0)| 

1/3- 



= 7r/3/9 in Eqns. (jSISI) . 
We observe that the approximate expression strongly de- 
pends on subleading power effects and is not a very good 
approximation for the exact result. Upon integration, 
these differences get reduced due to accidental cancella- 
tion in the intermediate /? region where the approximate 
results are smaller/larger than the exact one. 

Before closing this section we briefly explain our calcu- 
lational approach. The two-loop virtual corrections are 
taken from 0] , utilizing the analytical form for the poles 
Q. The one- loop squared amplitude is computed in 
The real-virtual corrections are derived by integrating 
the one-loop amplitude with a counter-term that regu- 
lates it in all singular limits [13] . The finite part of the 
one-loop amplitude is computed with a code used in the 
calculation of pp — > tt+jet at NLO ll|. The double real 
corrections are computed in [l2|. As in Ref. 12 1, we do 
not include the contribution from the reaction qq — >■ tiqq 
where the final state light pair has the same flavor as 
the initial state one. We expect the contribution from 
this reaction to be negligible. We consistently modify 
the collinear subtraction to account for this missing con- 
tribution. Details will be presented elsewhere. 



NUMERICAL PREDICTIONS: THE TEVATRON 

The NNLO results computed in this paper make it 
possible to predict the total top-pair cross-section at the 
Tevatron with high precision. Implementing the new 

13l | (with rrit — 



NNLO results in the program Top++ 
173.3 GeV and MSTW2008nnlo68cl pdf set [ij through- 
out) and adopting the scale and pdf variation procedures 



of Ref. Q we get (in pb): 



,NNLO _ 7 nnr: -^0.202(2.9%) r„„„, i +0.170(2.4%) r 

i .UUO n'3in/AA(r^\ [oCtlit&J _o 122 (1 7%) Lt^^ J ' 



-0.310(4.4%) 



fll) 



The fixed order NNLO result atJ^ includes the com- 
plete NNLO qq contribution ([3|) and the approximate 
NNLO result for the gg reaction as implemented in 



0|[l3|. We have set the unknown constant C, 



(2) _ 



0. 

We have verified that the sensitivity to the exact value 
of Cgg^ is around ±0.5% when Cgg' 



is varied in the range 
±1137, i.e. the value of the constant Cgg^^ [5|]. 

Our best prediction a^l = o-NNLO+nnll derived 
with full NNLL resumaton fis!] matched to the exact 

("2") _____ 

NNLO result for the qq reaction, including H^^' (|10l) : 



7.067 



hO. 143 (2.0%) 
-0.232 (3.3%) 



[scales] 



+0.186 (2.6%) 
-0.122 (1.7%) 



:pdf] , (12) 



while the contribution from the gg reaction to Eq. (fl^ 
is implemented as NNLOapprox+NNLL [l[ll. We take 
v4 = 2 for the value of the constant A |16l| that controls 
power suppressed effects. 



We find a 0.4% sensitivity of cr, 



NNLO+NNLL 



to the value 



of the constant A. To be conservative in our error esti- 
mate, we exclude the scale dependent term at the level 
of the unknown two-loop constant in the gg reaction, see 
[isj . Their inclusion lowers the scale uncertainty from 
±2.7% to ±1.7%. On the other side, including these 
scale dependent terms brings about 0(1%) sensitivity 
to the value of the unknown hard matching coefficients 
^gg,i ' ^gg,s which offsets the reduction in scale variation. 

We conclude that the error estimate of our best result 
Eq. (fT^ takes into account all dominant sources of the- 
oretical uncertainty, and that the missing NNLO contri- 
butions from other reactions will affect the above results 
at the percent level, i.e. they are accounted for by our 
theoretical uncertainty. 

On Fig. [3] we present the dependence of our best pre- 
diction cr^NNLO+NNLL ^^le value of the top mass. It in- 
cludes scale and pdf variation. We find very good agree- 
ment with the latest measurements from the Tevatron 
ITj and note that the total theoretical uncertainty is only 
about one-half of the total experimental one. 



CONCLUSIONS AND OUTLOOK 

In this work we calculate the genuine NNLO correc- 
tions to the total inclusive cross-section for qq —> tt + X. 
After extracting the two-loop hard matching constant 
from our result, we augment the NNLO evaluation with 
soft gluon resummation with full NNLL accuracy. As an- 
ticipated, our NNLO+NNLL numerical prediction for the 
Tevatron has substantially improved precision in com- 
parison with NLO+NNLL or approximate NNLO results. 
Most importantly, the accuracy of our theoretical predic- 
tion exceeds the accuracy of the best currently available 
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FIG. 3: Dependence of the NNLO+NNLL cross-section 
^NNLO+NNLL^ Eq. OH the pole mass of the top quark: 

scale variation (white band); scale + pdf variation (red band). 



FIG. 4: Scale variation at LO (orange), NLO (blue), NNLO 
(red) and resumed NNLO+NNLL (black hnes). 



measurements from the Tevatron. We are confident that 
our results will provide new insight to the forthcoming 
Tevatron analyses at full dataset, and will help scruti- 
nize the Standard Model to a new level. 

The very high precision of our result will allow critical 
comparisons between different pdf sets as well as extrac- 
tion of the top quark mass with improved precision. It 
is also a step in the derivation of the dominant missing 
SM corrections to ApB, whose calculation through order 
0{ag) will be the subject of a forthcoming publication. 

In a broader context, given the small number of ob- 
servables known at NNLO, it is interesting to address the 
question of the convergence of the perturbative series for 
this observable. In Fig. (|4]) we plot the scale variations of 
the LO, NLO, NNLO and NNLO-I-NNLL approximations 
as functions of the top mass. Each approximation is cal- 
culated with a pdf of corresponding accuracy. We observe 
a significant and consistent decrease in the scale depen- 
dence with each successive approximation. The overlap 
between the scale bands of the successive approximations 
also indicates that our scale variation procedure performs 
consistently well at each perturbative order. 
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code for the evaluation of the one-loop virtual corrections 
in qq — ^ ttg [llj. and Z. Merebashvili for clarifications re- 
garding Ref. [9|. The work of M.C. was supported by the 
Heisenberg and by the Gottfried Wilhelm Leibniz pro- 
grammes of the Deutsche Forschungsgemeinschaft, and 
by the DFG Sonderforschungsbereich/Transregio 9 Com- 
putergestiitzte Theoretische Teilchenphysik. 
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